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Abstract 

A precise measurement of the angle a in the CKM triangle is very important 
for a complete test of Standard Model. A theoretically clean method to extract 
a is provided by —> pn decays. Monte Carlo simulations to obtain the BTeV 
reconstruction efficiency and to estimate the signal to background ratio for these 
decays were performed. Finally the time-dependent Dalitz plot analysis, using the 
isospin amplitude formalism for tree and penguin contributions, was carried out. It 
was shown that in one year of data taking BTeV could achieve an accuracy on a 
better than 5°. 


^ untimely deceased 

^ corresponding author, email: Alexander.Vasiliev@ihep.ru 


Preorint submitted to Elsevier Science 


4 February 2008 





1 Introduction 


The Standard Model (SM) [1] which incorporates the quark mixing Cabbibo-Kabayashi- 
Maskawa (CKM) mechanism [2] has been increasingly successful, supported with many 
precise experimental results. This strongly indicates that, at low energies, the SM is the 
effective description of Nature. However, there are reasons to believe that there exists 
physics beyond the SM. For example, from the astrophysical point of view, it is a serious 
problem that the matter-antimatter asymmetry in the Universe can not be explained 
solely from the CP violation in the SM, which originates from quark flavor mixing. This 
observation, together with others, leads one to believe that there is a new physics, most 
likely, at the TeV energy scale. One of the critical measurements to test the SM or to 
obtain strong indications of new physics are precisie measurements of the angles in the 
Unitary Triangles (UT), which are of non-zero area if CP violation exists. 

A unique program that would have allowed one to challenge the SM explanation of CP 
violation, mixing and rare decays in the b and c quark system was proposed by the BTeV 
project [4] at the Tevatron at Fermilab. The design of BTeV exceled in several crucial 
areas including: triggering on decays with purely hadronic hnal states, charged particle 
identihcation, excellent electromagnetic calorimetry and excellent proper time resolution. 
Exploiting the large number of 6 ’s and c’s produced at the Tevatron collider, the exper¬ 
iment would have provided precise measurements of SM parameters and an exhaustive 
search for physics beyond the SM. The complete physics objectives of BTeV included 
measuring the CP violating angles a, (3, and 7 of the UT. In particular, the measurement 
of a is difficult due to small overall rates and because the gluonic penguin rates are of the 
same order as the tree rates, causing well known difficulties in extracting the weak phase 
angle. Quinn and Snyder [5] have suggested a theoretically clean method to extract a 
from decays of the type pn. The hnal state of these decays is not a CP eigenstate, 

which results in the need of a Dalitz plot analysis. We focus on the measurement of 
a, via collecting a large sample of decays. Direct measurements from the 

B-factories demonstrate that the average hh decays (where h = (p, tt)) is known to a 
precision of 0(10°) [3], with the use of isospin. In this paper we demonstrate that BTeV 
could have done the measurement with a much better precision. 

The paper is organized as follows. Section 2 provides brief introduction of the CKM matrix. 
In section 3 we give a general overview of the BTeV project. In sections 4 and 5, respec¬ 
tively, we report on the expected reconstruction efficiencies and the Signal/Background 
ratio of the B^ (p^r)'^ decays in BTeV. Section 6 covers the phenomenological formalism 
of B^ (p7J')° decays. In section 7 we describe results of the time-dependent Dalitz plot 
analysis of the simulated B^ (p^r)® decays in BTeV. 


2 The CKM Matrix and the Angle a 


In the SM there are three generations of leptons and quarks. The physical point-like 
particles that have both strong and electroweak interactions, the quark, are mixtures 
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of weak eigenstates, described by a 3x3 unitary matrix, called the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix [2], 


/ 



htid ^us ^ub 
Vcd Vcs Vcb 





\ Vtd Vts Vt, ) 



( 1 ) 


The unprimed states are the mass eigenstates, while the primed states denote the weak 
eigenstates. The Id/s are complex numbers that can be expressed in terms of four inde¬ 
pendent real quantities. These numbers are fundamental constants of Nature that need to 
be determined from experiment. In the Wolfenstein approximation the matrix is written 
as [ 6 ] 


VcKM — 


1 - AV 2 
-A 1 
AA ^(1 — p — irj) 


-AX'^ 1 


( 2 ) 


where (A, A, p, p) are four mixing parameters with A = \Vus\ ~ 0 . 22 , A k, 0.8 (measured 
using semileptonic s and b decays [7]), and p represents the CP violating phase. This 
expression is accurate to order A^ in the real part and A^ in the imaginary part. It is 
necessary to express the matrix to this order to have a complete formulation of the physics 
we wish to pursue. 

The unitarity of the CKM matrix leads to various relations among the matrix elements 


KdC + = 0 

(3) 

VusV:, + + VtsV:, = 0 

(4) 

VubV:, + 14 ,k;, + Vt,v:, = 0 

(5) 


that can be geometrically represented in the complex plane as triangles. These are “unitary 
triangles”, though the term “unitary triangle” is usually reserved only for the bd triangle 
in (5) where the angles are all thought to be relatively large. This CKM triangle is depicted 
in Figure 1. It shows the angles a, (3, and 7 . These angles are dehned as 


a = arg 


VtdVtl 

VudV:, 


(3 = arg 


VcdV:, 

VtdV,l 


7 = arg 


VudV:, 

VcdV:, 


and can be determined by measuring CP violation in B decays. 


They can roughly be divided in two classes : 
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Fig. 1. Unitary Triangle. 

• decays that are expected to have relatively small direct CP violation and hence are 
particularly interesting for extracting CKM parameters from interference of decays 
with and without mixing. 

• decays in which direct CP violation could be signihcant and therefore that can not be 
cleanly interpreted in terms of CKM phases. 

B decays used to extract f3 belong to the hrst group, whereas decays that have been 
considered to measure a belong to the second one. 

The primary source for measurements of sin(2/5) are the decays of the type b —> ccs. The 
most statistically signihcant measurements of CP violation in the B system were made by 
BABAR [9] and BELLE [10], resulting in the average value of sin(2/5)=0.725±0.037 [11]. 

Measuring a is more difficult than measuring (3 in several respects. First of all, the decay 
amplitudes are modulated by Vub rather than Vcb making the overall rates small, of the 
order of 10“^ to 10“®. Secondly, the gluonic penguin rates are of the same order as 
the trees causing large theoretical uncertainties in cleanly extracting a from asymmetry 
measurements, alone. 

The decay B^ tt+tt” has been proposed as a way to measure sin(2Q;). However, the 
penguin pollution is quite large and cannot be ignored. Gronau and London [12] have 
shown that an isospin analysis using the additional decays B~ 7r“7r° and B^ 
can be used to extract a [13], but the hnal state is extremely difficult to detect 
in any existing or proposed experiment. B tttt has been seen but there is no B 
decay vertex information, therefore there is no way to perform a time-dependent CP 
violation measurement. If fact, the data that does exist has been used to limit the Penguin 
contribution to these decays, but the limit is not very restrictive. Lipkin, Nir, Quinn 
and Snyder [14] have extended the analysis in [12] to include other decays, among them 
B — pn. Snyder and Quinn [5] subsequently extended that work and proposed not only 
an isospin analysis, but a full, time dependent, Dalitz plot study of B ^ pn decay to 
measure a. 

1 

A sample Dalitz plot is shown in Fig. 2. A striking feature of this Dalitz plot is that the 
events are concentrated close to the kinematic boundary, especially in the corners. This 
kind of distribution is good for maximizing the interference, which helps minimize the 
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errors. Furthermore, little information is lost by excluding the Dalitz plot interior, a good 
way to reduce backgrounds. 



Fig. 2. The Dalitz plot for i?° pir ^ 7r~^TT Tr°. 

Snyder and Quinn have performed an idealized Dalitz plot analysis that uses 1000 or 
2000 flavor tagged background free events. Trials using 1000 events usually yield good 
results for a, but sometimes do not resolve the ambiguity. With the 2000 event samples, 
however,the ambiguities disappear. 

Recently BABAR has made an important step to improve the constraints on a, via 
studying p'^p~ decays [15]. Using the isospin analysis they determined that 

the solution compatible with the Standard Model is a = (100 ± 13)". The estimate is 
based on 232 millions T(4S') —BB decays. This mode has potential show stoppers to 
improving errors on a. The analysis assumed a 100% longitudinal polarization of the 
B^ —i> p'^p~-, if this is not true, an angular analysis is needed and requires a lot more 
data. However, BABAR measured the longitudinal polarization fraction fi = 0.978 ± 
0.014(stat)]'^o;o29(-s|/st) [15] which is consistent with one. 

Regardless of that the p7i system remains theoretically the cleanest way to extract a. 
Recently BABAR has performed first full time dependent Dalitz plot analysis [15] and 
extracted a = (113]'^i7 6)". In the following sections we will show how, with the BTeV 

detector, one would have signihcantly improved sensitivity on a using the full time- 
dependent Dalitz plot analysis of the B ^ pn decays. 


3 The BTeV Concept 


BTeV was designed as a second generation experiment to study CP violation in B decays. 
It would have made possible to carry out practically all measurements of CP violation and 
decays of the B hadrons accessible the asymmetric R-factories and at CDF and DO running 
at Tevatron and it could have done those measurements at a much higher precisions. The 
detector design is ideally suited to study B decays containing neutral particles, especially 
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the modes of interest here, (p7r)°. Towards the end of the decade, LHCb will go into 

operation with similar capabilities for all-charged states, but will not have a high-quality 
calorimeter or a trigger which is efficient for all B decay modes. 

The studies presented in [16-20] indicate that the forward direction at the Tevatron 
presents a number of striking advantages. First of all, there is a large cross-section for 
the production of correlated bb pairs. Secondly, the B hadrons that are formed have 
relatively large momenta, on average 30 GeV/c, and, therefore, their decay products 
do not suffer much from multiple Coulomb scattering. This would allow BTeV to make 
precision measurements of the spatial origins of particles and as result BTeV would be 
able to determine if they arose from B hadrons that traveled on the order of several mm 
prior to their decays. Furthermore the geometry was very natural for certain aspects of 
detector technology that signihcantly enhance the physics performance. For these reasons, 
the BTeV collaboration designed a detector with “forward coverage.” 

The physics case for BTeV involves reconstructing a variety of different decay modes of 
the B, Bg, and other b hadrons and, in many cases, following their time evolution and 
tagging the flavor of the parent B at production and at the moment of decay. These decay 
modes may involve charged hadrons, charged leptons, photons (prompt or from 7r°’s), and 
tertiary vertices from the b ^ c decay chain. The product branching fractions of many 
decay modes of interest, including any tertiary decays, are quite small, typically 10“® to 
10“^. This, together with the large background of minimum bias events, demanded that 
BTeV be able to reconstruct multibody hnal states, with a good resolution in invariant 
mass, and to handle very high data rates. In order to carry out the physics program, 
the detector must have the ability to separate decay vertices from the primary interaction 
vertex and to reconstruct secondary B vertices and daughter charm vertices. This requires 
a precision vertex detector. It must also be able to measure the time evolution of decays 
for time-dependent asymmetry studies. The most demanding requirement is to be able 
to follow the very rapid oscillations of the Bg meson in order to study mixing and CP 
violation. The detector must have the ability to distinguish pions, kaons, and protons 
from each other to reduce confusion among similar decays such as i? —> tttt and B ^ Kn 
so that decays of interest will not be contaminated by other decays, causing the resulting 
measurements to be diluted. Many key decay modes have 7r°’s, y’s, or particles that decay 
into them, such as p’s or p’s. Leptons, muons and electrons (positrons), appear in many 
key hnal states so good lepton identihcation is also required. Finally, many of the detector 
properties which are needed to isolate and reconstruct signals are also needed to perform 
“havor tagging.” 

The BTeV detector is shown schematically in Fig. 3. The covered angular region is from 
approximately 10 mr to 300 mr with respect to the anti-proton beam. When a B decay of 
interest is contained within the acceptance of the detector, there is a high probability that 
the decay products of the co-produced B will also be within the acceptance of the detector. 
Furthermore, since the charged B decay products are not degraded by multiple scattering 
in the detector material, which allows accurate determinations of B decay vertices. 

The key design features of BTeV include: 

• A dipole centered at the interaction region placing a magnetic held on the vertex de- 
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BTeV Detector Layout 



meters 


Ring Imaging 



Fig. 3. Layout of the BTeV Detector 

tector, allowing the use of momentum determination in the trigger. There are two open 
ends of the magnet. One open end allows particles to flow into the instrumented “arm.” 
The held is used by the tracking system to provide precise momentum determinations 
of all of the charged particles. 

• A precision vertex detector based on planar pixel arrays. The outputs are used in the 
trigger processor to hnd detached heavy quark decay vertices in the hrst level trigger. 
They also provide precise and unambiguous three-dimensional space points to help 
reconstruct charged particles; 

• Precision tracking using a combination of straw tubes and silicon microstrip detectors, 
inside the straws close to the beam line, where the charged particle occupancies are the 
largest. This system, when coupled with the pixels, provides excellent momentum and 
mass resolution out to 300 mr; 

• Excellent charged particle identihcation using a Ring Imaging Cherenkov Detector 
(RICH). The RICH provides hadron identihcation from 3-70 GeV and lepton iden¬ 
tihcation from 3-20 GeV, out to the full aperture of 300 mr, which is crucial since the 
muon detector and calorimeter do not cover the full solid angle covered by the RICH. 
The RICH has two independent systems sharing the same space. One has a gas (C 4 F 8 O) 
radiator and a Multianode Photomultiplier photon detector and the other has a liquid 
C 5 F 12 radiator and a Phototube photon detector. Both systems work in the region of 
visible light; 

• A high quality PbW 04 electromagnetic calorimeter with excellent energy resolution, 
position resolution and segmentation, covering up to 200 mr, capable of reconstructing 
hnal states with single photons, 7r°’s, ?7’s or and identifying electrons; 

• Excellent identihcation of muons out to 200 mr using a dedicated detector consisting 
of a steel toroid instrumented with proportional tubes. This system has the ability to 
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both identify single muons above momenta of about 10 GeV/c and supply a dimuon 
trigger; 

• A detached vertex trigger at Level 1 using the pixel detector information, which makes 
BTeV efficient for most hnal states, including purely hadronic modes. The trigger 
ignores low momentum tracks that have large multiple scattering and would thereby 
avoid creating false secondary vertices; and 

• A very high speed and high throughput data acquisition system which eliminates the 
need to tune the experiment to specihc hnal states. 


4 Reconstruction Efficiencies for decays in BTeV 


As stressed earlier, measuring the time dependent CP violating effects in the decays 
> (pvr)° —provides a theoretically clean way to determine the angle a of the 
Unitary Triangle, as shown by Snyder and Quinn [5]. We report on the expected perfor¬ 
mance of the BTeV detector for these decays, taking into account excellent reconstruction 
efficiency of the tt'^’s that is made possible with the electromagnetic calorimeter based on 
PWO crystals. 

Excellent mass resolution in the reconstruction reduces the background significantly, 
particularly near the edges of the Dalitz plot where the fm events lay. In addition, good 
resolution in the proper decay time is crucial to determine the angle a. 

The reconstruction efficiencies for B ^ pn were studied using GEANT3-based simu¬ 
lation [22]. We generated two samples: 250,000 and 250,000 p° 7 r°. Both samples 

were generated with a mean of two Poisson distributed non-beauty interactions per beam 
crossing. This number of interactions per beam crossing corresponds to running at the 
designed Tevatron luminosity 2 x 10^^cm“^s“^ and 132 ns bunch spacing. It should be 
mentioned that we used two separate Monte Garlo samples only to refine the selection 
procedure and determine the reconstruction efficiency and signal-to-background ratios. 
For the Dalitz plot analysis the interference between charged and neutral p-mesons was 
simulated. 

The analysis relies on BTeV event reconstruction software packages, including track 
reconstruction based on the Kalman filter method, vertex reconstruction, and shower 
reconstruction. 

With the use of the electromagnetic calorimeter, we would find many good 7 r° candidates. 
Photon candidates are required to have minimum reconstructed energy of 1 GeV and pass 
a shower shape cut designed to reject hadronic showers. We reduce the background rate 
by insuring that the photon candidates are not too close to the projection of any charged 
tracks to the calorimeter. 

Figure 4a shows a 77 invariant mass distribution of the B ^ pn events when the pairs 
have energy sum greater than 5 GeV and the vector sum of transverse momenta greater 
than 0.75 GeV/c^. The vr® signal is very clear; the 71 ^ mass resolution in this sample is 
3.7 MeV/cT 




Fig. 4. The vr^ signal in the 77 spectrum from —> p+vr” simulated events (a) and the efficiency 
of reconstructing the tt^ as a function of distance from the beam line (b). All events considered 
have at least one signal B decay. The sample includes a mean of two Poisson distributed 
non-beauty interactions per beam crossing. 


Candidate are two-photon combinations with invariant masses between 125 and 
145 MeV/c^. The 7r° reconstrnction efficiency depends on the distance from the beam 
line and is presented in Fignre 4b; the vr^’s are taken from the —>■ events; this 

simulation was run with a calorimeter larger than that proposed so we could view the 
dependence on radius. The denominator contains all events in which the event passes 
the trigger, all charged tracks are reconstructed, and the combination of the two charged 
tracks passes some vertexing and detachment cuts. In the calculation of the efficiency we 
use the “right” two charged tracks. However, if a pair of good charged tracks is combined 
with a background 7r°, it does not signihcantly increase the efficiency. 


We look for events containing a secondary vertex formed by two oppositely charged tracks. 
One of the most important selection requirements for discriminating the signal from the 
background is that the events have well measured primary and secondary vertices. We 
demand that the primary and the secondary vertices be well dehned by requiring dof < 
2 for their vertex hts. Once the primary and the B decay vertices are determined, the 
distance L between the vertices and its error ctl are computed. The quantity L/ai is a 
measure of the signihcance of detachement between the primary and secondary vertices. 
We require T/ctl >4. The two vertices must also be separated from each other in the plane 
transverse to the beam. We dehne rtransverse in terms of the primary interaction vertex 
position {xp,yp, zp) and the secondary vertex position {xs,ys, zs), namely rtransverse = 
{xp — xsy + {yp — ysY and reject events when rtransverse < 0.132 mm Finaly, to 
insure that the charged tracks do not originate from the primary, we require that both 
the 7 T~^ and the vr" candidates have an impact parameter with respect to the primary 
vertex (DCA) greater than 100 ym. 
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M(7l^7l-7lO),GeV/c2 

Fig. 5. Invariant mass 7r'*'7r“7r® (after cuts) for the simulated (a) and (b) p'^7r~ 

decays. Each event includes on the average two Poisson distributed non-beauty interactions per 
beam crossing, mixed to the beauty production interaction. The number of signal event were 
counted as the number of entries in the ±2cj interval around the B mass, minus estimated 
number of background entries. 



4 4.5 5 5.5 6 6.5 7 


M(7l^7l-7lO),GeV/c2 

Fig. 6. Invariant mass 7r'’“7r“7r® (after cuts) matched vs generator tracks and photons, for the 
simulated B^ p'^tt~ at 2 and 6 interactions per crossing. 



10 










Table 1 

Selection Criteria. The notation is defined in the text. 


Criteria 

Value 

Primary vertex criteria 

X' <2 

Secondary vertex criteria 

X'<2 

'^transverse (cm) 

0.0132 

Normalized distance Ljo 

> 4 

Distance L, cm 

< 5 

DCA of track, pm 

> 100 

tproper/^O 

< 5.5 

E^+, GeV 

> 4 

E^-, GeV 

> 4 

Pt(7r+), GeV/c 

> 0.4 

GeV/c 

> 0.4 

Isolation for 7 , cm 

> 5.4 

E^o, GeV 

> 5 

pt(7r°), GeV/c 

> 0.75 


< 0.06 

m.y^, MeV/c^ 

125 - 145 

GeV/c^ 

0.55 - 1.1 


When we calculate the invariant masses of the 7r’''7r“ and pairs, we require at least 
one of them to be compatible with the p mass, that is, between 0.55 and 1.1 GeV/c^. In 
addition, we use several kinematic cuts which reduce the background io B ^ pii without 
signihcantly decreasing the reconstruction efficiency. We require that p™™ divided by the 
scalar sum of the pt values of all three particles, be small. The vector sum 

psum jg fiegj^ed with respect to the B direction of flight which is calculated from the 
reconstructed primary and secondary vertices. We also make a cut on the B proper time 
decay requiring it to be less than 5.5 times the 5° lifetime {tproper/to < 5.5). The selection 
criteria are summarized in Table 1. 

The results are shown in Figure 5 for B^ p°7r° (a) and B^ p'^tt~ (b) Monte Carlo sam¬ 
ples, respectively. The B^ mass resolution in these samples is in the range 32-42 MeV/c^. 
The signal interval is dehned as ±2(T around the B mass, minus estimated background. 
The reconstruction efficiency is (0.18 ± 0.02)% for B^ p°7r° and (0.22 ± 0.02)% for 

BO p±7r^. 

Similar simulation studies were repeated with six non-beauty interations per crossing 
mixed to the beauty production interaction to estimate reconstruction efficiency for the 
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—> p+TT decay. Results were compared with those at two interactions per crossing. 
The statistics used to compare the two cases were 100,000 events. 

At six interations per beam crossing the mass resolution remains practicaly unchanged, 
as it is fonnd to be (44±3) MeV/c^. The B^ p'^7i~ reconstrnction efficiency is estimated 
at (0.2 ± 0.02)%. This represents the effect of only 10% as compared to two backgronnd 
interactions per crossing. 

However, the number of false Stt combinations that would pass the cuts appears to increase 
somewhat as the nnmber of non-beanty interactions per beam crossing goes np. To prove 
that most of the entries in the B^ mass region are trne 7 r~^7r~7r^ combinations coming 
from the B^ decay, we have done a check against generator level information. Results of 
comparison are presented in Figure 6. It is clear that the B^ signal dominates in both 
distributions; false Stt combinations could, in principal, mimic the signal but most of the 
Stt combinations are the correct ones. 

Using the previously calculated reconstruction efficiency we could expect to have ~ 1000 
flavor-tagged events and ~ 150 flavor-tagged p°7r° events per year (2 x 10^ s given 

that BTeV was assumed to run 10 months per year). The samples would include both B^ 
and B decays, with proper time measurements for both states. 

In principal, one can use the untagged sample in the likelihood (see Section 7) to extract a. 
Actually, this sample does not carry any information on a bnt it allows to extract other 
parameters related to direct CP violation and helps the fit converge. This leads to an 
improved resolntion on a when the nntagged sample is utilized. However, in this paper 
we present results obtained only with the tagged sample. 


5 Signal-to-Background Ratio in 5° —(p7r)° decays in BTeV 


The analysis by Snyder and Qninn [5] showed that with 2,000 backgronnd free events 
they could always find a solution for a. BTeV could have collected such a statistics within 
4x10^ seconds (approximately 2 years). But we expect some background whose effects 
need to be estimated. 

For a channel with a branching ratio on the order of 10“^ and efficiencies lower than 1%, 
it is necessary to generate at least 10^ bb backgronnd events. For this stndy we generated 
2 X 10^ generic bb events {B pn channels excluded) and processed them throngh 
the GEANT3-based full detector simulation. Each event contains a mean of two Poisson 
distribnted non-beanty interactions. Selection criteria listed in Table 1 are applied. To 
get the backgronnd estimate, we count all of the events between 5 and 7 GeV/c^, then 
we scale that nnmber down by the ratio of the signal region divided by the background 
selection region. 

The results of the analysis are presented in Figure 7a and b. The signal-to-background 
levels are approximetely 4:1 and 1:3 for and p°7r°, respectively, when there are on 

the average two interactions per crossing. 
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M(7i+7c-7uO),GeV/c2 (after cuts) 


Fig. 7. Background tt^ invariant mass for p^-K^ and p'^jr at 2 interactions 

per crossing (a,b), and B^ p^n~ at 6 interactions per crossing (c) 

We have also investigated the effect of a larger number of interactions per crossing 
on the p’''7r“ background, similar to the study on the signal sample. We merged our 
background sample with an additional sample of non-beauty eventsgenerated with a 
Poisson distributed average of four interactions per crossing. Charged tracks in the merged 
events were projected onto the calorimeter, and photons from both samples were added 
in. Thus, the full confusion of six interactions per crossing is simulated in the calorimeter. 
The way we did this study, the confusion is not present in the simulation of the tracking 
system but separate studies show that the charged particle tracking system is reasonably 
robust against six interactions per crossing. The analysis then proceeded as before. We 
have reprocessed 1.33 x 10^ events and have compared the results corresponding to this 
statistics at two or six interactions per crossing. The background to p'^7r~ at six 

interactions per crossing is presented in Figure 7c). We have found that at six interactions 
per crossing the background to p’''7r“ increased to 109 events, at compared to 56 

events at two interactions per crossing. This demonstrates that the background could 
increase at six interactions per crossing but the effect is expected to be about a factor of 2. 


6 Representation of amplitudes and phenomenological inputs 

6 .1 Classification of the amplitudes 

In this section we hrst dehne the formalism. Next, we estimate parameter values so that 
the simulation is as close as possible to reality. 

Amplitudes of neutral B^ meson decay to pn are represented in the form 

|R°) = fia,„ {tj} = {+-}, {-+}, {00}, (7) 
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( 8 ) 


Ctij 


(-e-'“ry + Py 


e 


-i /9 


1 


where Tij and Pij give tree and penguin amplitudes, correspondingly, as depicted in Fig. 
8 extracted from [23]. 



Fig. 8. The tree (left) and penguin (right) diagrams for the p+vr - ( 5 ° ^ p- ■7r+) decays. 


The fk represents the relativistic Breit-Wigner form for p tttt 

_ cos{0k) 

^ s — w?p + fn(s) ’ 


(9) 


where s is the square of the invariant mass ( 7 ri, 7 r 2 ) and 6 k is the angle between a decay 
pion and the line of flight of the p. The function n(s) is dehned as 


n(9) 


^ / p{s) 
v/i \p{ml) 


rK), 


p(s) = 


( 10 ) 


where nip is the p mass and F is the width. 

The amplitudes aij for and decay are written as a sum of tree (T) and penguin 
(P) contributions as 


a+_ = -P^T+- + e-*^P+- 
a_+ = -e*^T-+ + e-*^p-+ 

aoo = 

a+_ = -e-^^T-+ + 
a_+ = -e-^^T+- + pf^P+- 

aoo = _e-*9T“ + e*^P“, (11) 

where 7 and (3 are the usual CKM angles and q; + /5 + 7 = tt. Using both isospin symmetry 
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and the fact that the Penguin amplitude is a pure AJ = 1/2 transition leads to the 
replacement 


p«, = -j {p+- + p- 


( 12 ) 


For the tree diagrams after Fierz transformation one naively gets 

Too = i - (r+_ + T_+), 

2 ai 


(13) 


where factors ai ,2 represent contributions due to gluon corrections to the weak interac¬ 
tions of quarks and depend on renormalization scale. The corresponding Lagrangian is 
approximated by 


C = 


2^2 


VudV:,C^{h0^u^){ukO>^d 




± <5*5) 


(14) 


where indices mark SU(3)-colors, and the factors are dehned as 

Oi = ^ + 1) + C.(N, - 1)1, 02 = ^ lC+(Af, + 1) - CL(N, - 1)] (15) 

where = 3. In the limit of neglecting the gluon corrections we get 

I? ^1 I? ^2 /Y~* 

However, the corrections taking into account the renormalization group dependence on a 
decay scale put the ratio 02/01 to a negative value approximately given by 

— =-0.25 ±0.05. (16) 

Oi 


In fact, as was shown in [24] and [25], expressions (15) following from the factorization 
hypothesis can be signihcantly modihed by “non-factorizable effects” in the complex 
phase. This means that only the absolute value of the ratio 


02 

Oi 


0.25 ±0.05 


is reliable. The value of the phase is not realiably predicted by the theory but is roughly 
estimated to be about 45°. 

Corrections in (12) and (13) due to isospin symmetry breaking are considered to be 
negligible [23] and they are not included in the following simulations. 



6.2 Phenomenological constraints on the parameters 


A phenomenological analysis of possible values for the amplitudes in (7) has been per¬ 
formed in [23] and [26]. It is based on a global £t to the measured rates assuming 
SU(3)-flavor symmetry for B pn, B K*n and B pK decays. The analysis 
gives approximately twice enhancement of penguin amplitudes in comparison with QCD 
expectations [25]. The preferable values of amplitudes with theoretical expectations of 
uncertainties are summarized in Table 2, which has two sets of parameters we use in our 
modelling of the signal. 


Table 2 

Amplitudes in units of _set to 1 and expected uncertainties from fits in [23] and [26]. 


Parameter 


\T-+\ 

arg [r_+ 

arg[P+_ 

\P-^ 


arg 


T-+ 

'P-+ 


arg 


r_ 

ai 

ai 


a 


Set I 

Set II 

Theor. uncertainties oi 

0.8 

0.8 

0.63-0.9 

-20° 

1 

to 

o 

o 

±10° 

0.18 

0.18 

±0.05 

CO 

o 

o 

O 

O 

CO 

±30° 

0.28 

0.28 

0.14-0.32 

00 

o 

o 

130° 

±60° 

0.25 

0.25 

0.18-0.32 

45° 

45° 

0 - 27r 

00 

00 

o 

100 ° 

80° - 110° 


The value of |T__|_| is ordinary fixed by the factorization hypothesis [23], so that it is equal 
to the ratio of decay constants f^/fp ~ 0.63, while the ratio fitted by SU(3) ansatz results 
in a greater value of about 0.7. Nevertheless, we fix this parameter to 0.8, reproducing 
the mean magnitude of branching fractions. The uncertainties of complex phases are not 
given explicitly in [23] and [26], but we expect them to be lower than 30° at fixed absolute 
values of penguin-to-tree ratios. The amplitude of B^ p°7r° is constructed in accordance 
with (12), (13) and (16). 

Now we make the transition to the estimate of parameter values in the simulation and to 
the comparison with the existing experimental results. 

The overall normalization is tuned to the experimental sum of CP-averaged branching 
ratios 

= ey + Bp = (24.0 ± 2.5) X 10-®, 
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where 


B+- = i {bib” ^ p+TT-] + B|B“ ^ p-TT+l} = (13.9 ± 2,2) x IQ-”, 

i {B[B” ^ p-T!*] + B|B“ ^ p*T!~]} = (10.1 ± 2.1) X 10“”. 

For example, taking the set I of parameters we get 

B[B^ p+TT"] = 16.5x10-®, B[B^ p-7r+] = IFlxlO”®, B[B^ pV] = O-OxlO”®, 
and 

^ p+TT"] = 4.6x10-®, i3[.B® ^ p-7r+] = 11.6x10-®, ^ p®7r°] = 1.4x10-®, 

which should be compared with experimental averages from BELLE, BABAR and CLEO 
in [23] and [26] 

B[B^ p+TT-] = (16.5t|^) X 10-®, B[B^ p-7r+] = (15.4l|2) x Ip-®, 

and 

B[B^ p+TT-] = (4.8l|^) X 10-®, B[B^ p-7r+] = (11.4+|®) x 10-®. 

The above branching ratios give CP-averaged values of 

B+- = 14.0 X 10-®, = 9.3 X 10-® 

for the charged modes, while for the neutral mode we have 

= 1 {i3[R0 ^ ^ _ ^_2 X 10-® 

are consistent with the experimental value 

i3®® = (1.7 ± 0.8) X 10-® < 2.5 X 10-® at 95% C.L. 

Note, that the neutral CP-averaged mode weakly depends on the complex phase of 02 / 01 , 
but branching ratios of R® and R® strongly depend on that phase: for instance, putting 
arg[o 2 /oi] = TT gives B[B^ p®7r®] = 2.1 x 10-® and B[B^ p®7r®] = 0.2 x 10-®. 

The time-dependent CP asymmetry is given by 
r(:B°(t) - p%^) - V{B\t) - p%^) 

CL/^p — _fi — 

T{B (t) p^TT^) + r(5^(t) ^ p^TT^) 

= {Sp^ ± ASp^) sin(Amdt) - {Cp^ ± ACp^) cos(Amrff). (17) 

In this formula Sp^, and Cp^, represent mixing-induced CP violation and flavor-dependent 
direct CP violation, respectively. The value of AS'p^ and ACp^ are CP-conserving. The 
ACpj, characterizes the asymmetry between rates r(i?® —p’''7r-) -|- r(i? — p~'n^) and 
r(R® —> p-vr"'') -|- r(i? —p’'"7r-) at t = 0, i.e. at initial moment of evolution, while AS'p^ 
indicates mixing of decays at t % 0, and as we have found, it strongly depends on both 
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the relative strong phase of penguin with respect to tree amplitude in p mode (the 
parameter arg[P_+/T_+]) and the CKM angle a. 

Time-integrated asymmetries are given by 


dip. + gp. + _ iV(P° ^ p-7r+) - 7 V(eo ^ p+TT-) 

1 + ACp. + dlp.C'p. iV(;B° ^ p-7r+) + iV(P0 ^ p+TT-) ’ 


dip. - Cp. - Ap.AC'p. N{W ^ p+TT-) - iV(po ^ p-7r+) 

1 ~ ACp. — Ap^Cp. N{b'^ + N{B^ 


|a+_| 

^ + 0+- 

2 _ 

a_ 

2 _ 

a_ 

2 

|a+_ 

P + h+_ 

P + 

a_ 

P + l 

a _ 

T 


(18) 

(19) 

( 20 ) 


where non-zero values of Ap^ and ApJ" indicate direct CP violation. 

In Table 3 we show a comparison between observable quantities obtained by using sets I 
and II values with experimental data. 

The comparison shows that the range of parameters we use seem to be reasonable. 


7 Time-dependent Dalitz Plot Analysis of > (p^r)® decays in BTeV 


In section 4 we demonstrated that BTeV would be able to collect a sample of ~ 1000 
flavor-tagged B^ — {pn)^ events within one year of operation, which would allow a reliable 
Dalitz plot analysis of this decay mode. 

The model of the Dalitz plot analysis has three parts : 

• B pTi signal 

• random true p plus random true tt, the “resonant background” 

• uniform density, the “non-resonant background” 

The formalism used to £t the Dalitz is based on 13 independent parameters : 6 amplitudes, 
6 strong phases and the weak phase a itself. Using the constrains given in eq. 11 and 12 
we can reduce the number of parameters to 11. We fix a reference rate and strong phase 
so that the total number of parameters reduces to 9. Two additional parameters must be 
added if we allow the resonant and non-resonant background fractions to be determined 
by the fit. 

Due to the low reconstruction efficiency of this particular final state it would not be 
feasible to include the Snyder-Quinn formalism directly into the full detector Monte Carlo 
simulation; it would have required significant computer power and the generation of a huge 
number of events to obtain the desired statistics. 



Table 3 

A comparison of values of quantities evaluated from sets I and II with available experimental 
data from [26]. 


Quantity, units 

Set I 

Set II 

Exp.data 

B[B^ - 

P+TT ], 10 ® 

16.5 

16.5 

16.5l|i 

B[B^ - 

p“7r+], 10“® 

14.1 

14.3 

15.01|2 

B[B^ - 

p°7r°], 10-® 

0.9 

0.6 

1.7 ±0.8 (^00) 

1 

O 

SQ 

P~^TT~], 10“® 

4.6 

6.6 

4.8l|6 

B[B^ - 

p“7r+], 10“® 

11.6 

11.4 

11.6l|8 

B[B^ - 

p°vr°], 10“® 

1.4 

1.8 

1.7±0.8(SOO) 


A 

^pTT 

-0.100 

-0.054 

-0.114 ±0.067 


Sp-n 

-0.15 

-0.30 

-0.13 ±0.18 


ASp^ 

0.33 

0.35 

0.33 ±0.18 


CP'K 

0.33 

0.27 

0.35 ±0.14 


ACp^ 

0.24 

0.16 

0.20 ±0.14 



-0.17 

-0.18 

-0.18 ±0.14 


A —^ 

-^pTT 

-0.51 

-0.37 

-0.52l°iS 


a 

88° 

100° 

ioo°ii2:, 





(CKM unitarity; 98° ± 16°) 


We have opted for a different approach. The generated template events are distributed flat 
over the Dalitz plot domain, with the exponential time distribution and random tag=±l. 
We further use a rejection algorithm based on the isospin amplitudes formalism for the 
tree and the penguin contributions to the Dalitz plot. 

Time evolution of the (p7r)° decay amplitudes, including 5 — 5-mixing, is given by : 


p./n. AMt r. .q . AMt 
A = e"^*/^(cos- 5°) -h t- sm- 

2 p 2 




A = e ' (^-sm-15 )-I-cos- \B )) 


( 21 ) 

( 22 ) 


where |5'^) is given by eq.7. The template events are accepted or rejected based on whether 
a random number is less than or greater than |^P/|v4.maa:P- 

The background has been parameterized to account for both non-resonant and resonant 



components. The non-resonant background has been uniformly distributed over the Dalitz 
plot domain. The resonant background allows the two pions to have a Breit-Wigner shaped 
enhancement with the p line shape. 

The process of reconstruction of the accepted events is simulated by smearing them using 
the resolutions on momentum reconstruction and lifetime. These values were obtained 
from the simulation described in Section 4. The smearing has been computed comparing 
the reconstructed momentum, of 7r+, tt~ and 7r° to the generator information, a{pgen — 
Prec)/Pgen '■ 0.7% for charged pions and at 0.9% for tt ^’ s . Signal events are generated with 
an exponential time distribution. The rejection algorithm appropriately shapes the time 
evolution of the according to mixing and CP violation. The resolution on lifetime has 
also been estimated to be 64 fs using the Monte Carlo described in section 4, by computing 
the reconstructed lifetime and comparing it to the generated one. It should be pointed 
out that the resolution on lifetime is independent of lifetime. Proper time dependent 
acceptance was included in the likelihood. The background level is determined by a full 
Geant simulation of 20,000,000 generic bb events; it is assumed that this background has 
an exponential time dependence given by the average lifetime of 6-flavored hadrons. 

We have used two values for a : a = 88° and a = 100°. For each case we have generated 
500 independent trials, starting with different random numbers every time. Every trial 
contains 1000 signal events, 250 non-resonant background events and 250 resonant back¬ 
ground events. This corresponds to one year (2 x lO'^s) of data taking. The background 
level was chosen based on the MonteCarlo studies described in section 5, for the case of 
of running at 396 ns bunch spacing (6 interactions per beam crossing), which would be 
BTeV’s most challenging scenario. The Dalitz plot for one such sample is shown in Fig. 9. 


To extract the parameters and the associated errors we have used an unbinned maximum 
likelihood £t. The likelihood over the full Dalitz domain is given by 


-2lnC = —2 /n 






2=1 


1 -|- Rn 


Rr 


2 /ra 
i=i 




1 T Rnon T Rres 


(23) 


where 




+ 


R x^ + R X 

^ ^non K r ' Wes ^ 

Rft Mbw 


X 


|vA(Sj , Sj , tj, CK, .. 




A/'(a,..) 


Rnon X Rres X 


X £(Sj', Sj ) + 


Ml 


BW 




(24) 


(25) 
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Fig. 9. Dalitz-plot and proper time distribution for one trial of 1000 events for a = 88° (detector 
efficiency included). 

Here = (mj,+ +m^o)^ and s~ = (rriT,- +mj,o)j are two Dalitz plot variables for the j-th 
event. N^o and N-^ are the total nnmber of the and events. The normalization 

M ={\A\‘^ + \A\^) X e is integrated over the Dalitz plane and over the proper time and 
weighted by the detector efficiency. Rnon and Rres are dehned as 


l^jsignal 


jSJback 

res 

■\Tsignal 


The £t has been performed for 500 samples, for each inpnt valne of a = 88° and a = 100°, 
to confirm the stability of the procednre. The resnlts of the fit valne for a (right) and 
its deviation from the inpnt valne (left) are presented in Fig. 10, for the inpnt valnes 































Table 4 

The change (in degrees) of the value of the a coming out of the likelihood fit at different errors 
in the tagging dilution factor (in %). 


input a, degrees 


2 

5 

10 

15 

20 

25 

88 

additional displacement, degrees 

+1.0 

+0.4 

-0.6 

-1.6 

-2.7 

-3.9 


-0.9 

-1.6 

-2.9 

-4.3 

-5.8 

-7.2 


of a = 88° (Sample 1, Fig. 10a) and a = 100° (Sample 2, Fig. 10b), respectively. We 
measured a = (88.3 ± 1.6)° and Aa = (1.7 ± 0.09)° for Sample 1 and a = (99.7 ± 1.5)° 
and Aa = (1.8±0.1)° for Sample 2, which is in good agreement with the input parameters. 

These results have been obtained for ideal tagging. We have also made hts for different 
errors in the tagging dilution factor. The results are presented in Table 4 in the range 
from almost ideal adu = 2% to the very conservative adu = 25%. The displacements can 
be considered as systematic errors and should be summed with the ideal tagging error in 
quadrature. Assuming that the error in the tagging dilution factor in BTeV would have 
been 10-15%, we estimate the accuracy in measuring a at (1.8 — 2.3)° for a = 88° and 
(3.4 — 4.7)° for a = 100°. We consider the most conservative case, a 15% error on the 
dilution factor for a = 100°, and conclude that BTeV could measure a with the accuracy 
of better that 5° in one year of operation. 

In order to ensure that our event generation model is correct (see eq. 17-20) we have 
integrated eq. 23 over the Dalitz-plot domain, so that only the proper time dependence of 
the 5-meson decay rate is left. We £t the data with the assumption of the time-dependent 
CP asymmetry of the 5-meson decay as expressed in eq. 16. The CP asymmetries a'^p for 
p’'‘7r“ and a“ for (see eq. 17) are shown in Fig. 11 : the solid curves represent the 

£t results. The results obtained for the parameters in eq. 17 are : 


Sp^ = -0.22±0.06(-0.15) 

(27) 

ASp^ = 0.29 ±0.06(0.33) 

(28) 

Cp^ = 0.38 ±0.05(0.34) 

(29) 

ACp^ = 0.24 ±0.05(0.23) 

(30) 


where the numbers in parentheses represent the input to the simulation. 

The results of the £t are in good agreement with the input values used in the Monte Carlo 
simulation. This justihes the validity of the model we used to extract a and demonstrates 
that these important parameters could have been measured with the accuracy of 0.05-0.06 
in one year of BTeV operation. As one can see in Table 3, up until now these values are 
known only crudely. 
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(b) Input value a = 100° 

Fig. 10. Accuracy in determining a : fit value of a and deviation from input value {ainput — ctfu) 

8 Conclusions 

Physics simulations of the decay B ^ p7i for the BTeV project at Fermilab has been 
performed. The main idea was to estimate the expected accuracy in extracting the angle 
a of the UT. 

To calculate the signal-to-background ratio for the decay of interest, 2 x 10^ background 
events were simulated and processed through the full detector simulation based on GeantS 
package. Using the excellent electromagnetic calorimeter based on Lead Tungstate crys¬ 
tals, the 5° —p’''7r“ decay the signal-to-background ratio is estimate at 4:1 or 2:1, for 
132 ns or 396 ns beam crossing intervals, respectively. 
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Fig. 11. CP-asymmetries for and a~p for The solid curves represent the fit. 


A phenomenological analysis has been made for the possible values of tree and penguin 
amplitudes and phases for the process of interest, based on a global £t with SU(3)-flavor 
asymmetry for B pvr, B K*7r and B pK. The latest experimental data from 
BaBar and BELLE were used in this analysis. 

Dalitz-plot analysis of the B ^ pn decay with input from the phenomenological analysis 
has been presented. It has been shown that in one year of data taking BTeV could achieve 
the accuracy better than 5° on the angle a. 

The interference between tree and penguin diagrams can be exploited by measuring the 
time dependent CP violating effects in the B ^ p7i decays. In this paper it has been 
found that mixing-induced CP violation parameter and flavor-dependent direct CP 
violating parameter could be measured with the accuracy of 0.05-0.06. 
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